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TECHNICAL NOTE 25^6 


VISCOSITIES OF AIR AND NITROGEN AT LOW PRESSURES 
By Herrick L. Johnston^ Robert W. Mattox, and Robert W. Powers 


SUMMARY 


Information is presented of the results of an investigation to 
measure the true viscosities of nitrogen and air at 306 ° and 273 ° K at 
pressures from 500 to 0.0005 millimeter of mercury and at 19 ^° and 79° K 
at pressures down to 0.3 millimeter of mercury. 

In the region in which the thickness of the gas layer is many times 
the mean free path of the gas (the region of slip), the "slip coeffi- 
cient" - defined as the proportionality constant existing between the 
slip velocity and the velocity gradient - was found to show a linear 
dependence on the inverse pressure down to about 0. 02 millimeter of 
mercury. 

In the region in which the pressures are sufficiently low to permit 
the neglect of collisions between molecules, theory has indicated that 
the viscous drag should be directly proportional to the pressure and 
inversely proportional to the square root of the temperature. Results 
of the present paper ihdicate that this relationship was reached at 
pressures a little below 0.02 millimeter of merctiry. 


INTRODUCTION 


Ordinary gas kinetic theory leads to the conclusion that the vis- 
cosity of a gas at low and moderate pressures is Independent of its 
pressure. This conclusion has been confirmed experimentally for pres- 
sures between about 0.1 and a little more than 1 atmosphere. At higher 
pressures, there is an Increase in viscosity, that in the case of carbon 
dioxide (CO 2 ), for example, amounts to about 8 percent at atmospheres 

and to about 200 percent at 100 atmospheres. At the latter pressure, 
the viscosity of CO 2 is almost lineeir with its density. (See refer- 
ence 1, p. 150 . ) 

At pressures below about 0.1 atmosphere, the meastired viscosity 
decreases with decreeuslng pressure, gradually approaching a linear 
relationship. Down to pressures of about 0.001 atmosphere, this decrease 
is uBvially attributed to "slip" of gas layers over adjacent surfaces. 
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It is intei^reted as being due to a shai^) decrease in the velocity 
gradient, normal to tbe surface, of gas flow parallel to the surface. 

(See reference 1, p. 292. ) There is an apparent reduction of viscous 
drag on the sxirface due to slip. This effect appears as an apparent 
reduction -in the viscosity of the gas, since viscous drag always occurs 
on some surface in viscosimetry. It is probable, however, that the 
gas maintains its ordinary viscosity at some distance from the surface, 
in the sense that there is no reduction in the shearing force between 
InininRr layers of gas that axe moving parallel to each other with 
different velocities. 

The measured values of the viscosity in the intermediate pressure 
region decrease with pressure and depend upon the dimensions of the 
apparatus, particularly iqjon the distance between the disks or cylindeis 
in an oscillating- disk or cylinder viscosimeter. Such measurements, 
however, can be used to obtain a coefficient of slip ^ which can then 
be used to compute the viscous drag of an apparatus of different 
dlmens ions . 

The conception of slip is only applicable in the interpretation 
of viscous drag at low pressxrres when the thickness of the gas layer 
is many times the mean free path of the gas. Men the pressure becomes 
so low that collisions ih the gas phase are Infrequent compared with 
collisions with the sirrface, neither the ordinary concept of a shearing 
force between laminar gas layers nor the concept of surface slip is 
applicable. However a viscous drag still exists and can be explained 
in terms of an exchange of momentum between the gas molecules and the 
moving surface. Values computed from viscosimeter data in the region 

of very low pressures, \isually below 10”5 atmosphere, are termed "free- 
molecule viscosities." They are directly proportional to the pressure 
and are independent of the dimensions of the apparatus, provided there 
is a hi^ proportion of surface collisions. 

Measurements that established the foregoing concepts of viscous 
drag at low pressures were made by Kundt and Warburg (reference 2), by 
Timiriazeff (reference 3)^ by Hogg (reference 4), by Stacy (reference 5)^ 
by Van Dyke (reference 6 ), and by Blankenstein (reference J). However, 
the total literature on low-pressure viscosities is relatively sceinty. 

An investigation was therefore initiated at this Laboratory to determine 
the viscosities of air and of nitrogen at low pressures and low tempera- 
tures. Apparent viscosities have been measured at 306° K and at 273° K 
over the pressure range 0.0005 to 500 millimeters of mercury (Hg) and 
at 194 ° K and at 79 ° K down to pressures of 0.3 millimeter of Hg. The 
latter measurements coiild not be extended to lower pressures because of 
erratic drifts of the "zero point" of the suspension in the oscillatlng- 
disk viscosimeter. 
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This work was done hy The Ohio State University under the sponsor- 
ship and with the financial assistance of the National Advisory Committee 
for Aeronautics. 


APPARATIB AND PROCEDURE 


The gases used in this investigation were carefully purified. The 
air was freed of carhon dioxide and of moisture hy passing it throu^ 
successive towers of concentrated potassium hydroxide (KOH), concentrated 
sulfuric acid ^H 2 S 0 ]j^ , and stick potassium hydroxide (KOH), and thence 

throTogh a dry- ice trap and over a drier of phosphorous pentoxide ^P 20 ^^. 

Hy-Purity nitrogen was obtained from the Linde Air Products Co . , who 
represented it to he better than 99-99 percent pure, and was passed 
throu^ a dry- ice trap and over Measiirements were made in an 

osclllating-disk viscosimeter described in an earlier paper (reference 8 ) 
from this Laboratory. It has a moving disk of burnished silver, 1 milli- 
meter thick and 3 centimeters in diameter, that executes rotational 
oscillation between two fixed silver disks. There is a clearance of 
1 millimeter between the movable disk and each of the fixed disks. 

The constantan suspension wire used in the original apparatus 
(reference 8) was replaced by a 0,001- inch manganin wire, in order to 
reduce the energy dissipated in torsion of the wire. The method used 
to anneal constantan suspensions, heating electrically in vacuo to a 
dull red, was unsuccessful with manganin and resulted in a brittle 
manganin wire. The manganin wire was, therefore, annealed by heating 
at 150° C in a vacuum oven for 48 ho\irs. 

The pressure inside the apparatus was measured by two methods. A 
mercury meuiometer was used for pressures between 10 and 5 OO millimeters 
of Eg, and a McLeod gage for pressures below 10 millimeters of Eg. 

Experiments were conducted at 306° K by immersing the bulb of the 
apparatus in a manually controlled constant-temperature bath of water, 
at 273.2° K in a mixture of ice and water, at 195° K in a dry- ice - 
acetone bath, 6ind at 79° K in liquid nitrogen. 


CALIBRATION 


Maxwell (reference 9) derived the following fundamental equation for 
the oscillatingTdlsk viscosity apparatus; 

q = (X, - x)/Ct 


( 1 ) 
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where t] Is the viscosity^ X the logarithmic decrement^ x the wire 
constant, t the period of one oscillation in seconds, and C the 
apparatus constant. 

The logarithmic decrement can he represented as 



where Lq is the full amplitude of the nth oscillation and n takes 
on the values 1, 2, . . . , i. The present custom has been to take 
series of 10 oscillations (i = 5) for gas runs. The five values for X 
thus obtained ordinarily agreed to within 0.1 percent. 

As in the earlier work (reference 8), the apparatiis constant was 
determined hy measurements of logarithmic decrements in air at tempera- 
tures close to 296.1° K. The viscosity of air at 296.1° K was assumed 

to he 1833.0 X 10"7 poise. Wire constants were determined from log- 
arithmic decrements in vacuo. Values of the wire constants for the 
three different manganin suspensions used in the course of this investi- 
gation and the corresponding apparatus constants are given in table I. 


EXPERIMENTAL RISULTS 


The data obtained for nitrogen are given in table II and figure 1, 
and those for air in table III and in figure 2. The viscosities shown 
in the tables were con5)\ited by equation (l), using the values of the 
apparatus constant listed in table I. Corrections were applied for the 
changes in dimensions dvie to thermal contraction of the metal disks 
(reference 8) and, at low pressures and at temperatures below 273° K, 
for sli^t sijperheatlng of the gas between the disks due to heat leak 
down the suspension system. The latter correction amounted to about 
0.2 percent for air at 193° K, but Increased to about 1 percent at 80° K. 

Data obtained earlier at this Laboratory by Johnston and McCloskey 
(reference 8) were interpolated to the temperatures employed in this 
investigation. When the present results for pressures above 200 milli- 
meters of Hg are compared with those of reference 8, which correspond 
to atmospheric pressure, a mean deviation of -O.I7 percent exists for 
nitrogen and of 0.00 percent for air*. 
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DISCUSSION AND INTERPEETATION OF RESUIZTS 
Region of Slip 


Viscosity can "be defined as the coefficient of shearing force 
between laminar layers of fluid that move parallel to each other at 
different velocities. For layers in the IZ-plane that move in the 
Y-direction, this quantity can be represented by the relationship 

Py = ri (dVy/dX) (3) 

where Vy , is the mass velocity of the fluid and varies from layer to 
layer, dVy^dX is the velocity gradient normal to the direction of 
flow, and t] is the coefficient of the Y component of the shearing 
force P. 


In ordinary viscous flow parallel to a solid surface, there is 
a constant velocity gradient dV^-^dX normal to the surface. This 
gradient extends from the thin fluid layer adjacent to the surface, 
whose velocity relative to the surface is zero, to the point of the 
fluid where the velocity is at a maximum. This gradient produces a 
viscous drag on the surface. 

. With slip flow, the velocity gradient is uniform in the bulk of 
fluid but increases rapidly near the surface. The result is as thou^ 
the velocity itself decreased uniformly to a finite value u at the 
siirface and then fell abruptly to zero. The effect would correspond 
to slip of the bulk fluid over the surface at a velocity u. In terms 
of ordinary frictional forces, a proportionality should exist between 
the slip velocity u and the velocity gradient dVy^dX. This relation 

defines a proportionality constant termed the "slip coefficient," 
as 

u = ^(aVy^dx) (4) 

The consteint 5 depends upon the nature of the surface, particularly 
upon its degree of roughness. Probably it is also dependent upon the 
characteristics of the fluid. Equation (4) shows that ^ also corre- 
sponds to the distance by which the coordinate X must be extended 
beyond the surface to reduce Vy to zero, if the gradient remained 

constant at its slip value. Maxwell (reference 9) showed that the slip 
coefficient almost equals the mean free path, for slip over smooth 
surfaces . 
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It is possible to evaluate the slip coefficient from measurements 
with an oscillating-disk viscosimeter if the apparent viscosity, cal- 
culated from data taken in the slip region of pressure, is compared 
with the true viscosity obtained at normal pressures. This evaluation 
can be made by combining equation (l) with the following equation for 
the apparatus constant, by Keesom and Van Itterbeek (reference 10) ; 

c = + ai)/(>nau ^i) (5) 

where r is the radius of the rotating disk, I the moment of inertia 
of the suspension system, and d]^ and dy^ the spacings that separate 

the rotating disk from the upper and lower fixed disks, respectively. 
This equation reduces to 


C = itr^/2dl (6) 

for an apparatus like the present one, in which the rotating disk is 
centered between the two fixed disks so that dy = d-|^ = d. 

The apparatus constant In equation (l) is not a true constant when 
slip occurs. • Using eqmtion (l) with an apparatus constant determined 
in the viscous region, therefore, will not give the true viscosity but, 
instead, will give an apparent viscosity that decreases with pressures. 
HOTrever, the true gas viscosity, which is associated with the uniform 
velocity gradient away from the wall, can be calculated from eqmtion (l) 
if the apparatus constant is corrected for the reduction in viscous drag 
due to slip at the surfaces of the disks. This correction is made by 
adding 2^ to the plate separation d of equation (6). If T|p repre- 
sents the apparent viscosity obtained from the uncorrected form of 
equation (l) and if t] is the true viscosity of the gas obtained by 
substituting the corrected apparatus constant for slip Cg into equa- 
tion (l), the following relationship is obtained: 

ilAp = C/Cg = (d + 2^/d 
which can be rearranged as 

5 = (d/2)||Tl/Tlpj - 2 ^ (7) 

Equation (7) has been used to calculate the slip coefficients of 
air and of nitrogen from the present data for pressures between 10 and 
0.01 millimeter of Hg. The results are summarized in tables IV and V 
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and are plotted in figures 3 and 4, which show a linear dependence on 
the inverse pressirre. Constants in the equation 

t = K(l/P) (8) 

are listed in table VI for both air and nitrogen. Figure 5 is a plot 
for nitrogen at 306° and at 273° K extended to 0. 01-millimeter pres- 
sure. As the greph indicates, the linear dependence on l/P continues 
down to about 0.02 millimeter of Eg where noticeable departure from 
linearity begins. 

It was thou^t, from the relationships shown in table VI, that 
the slip coefficients mi^t be proportional to the specific volumes 
or reciprocals of density at all temperatures. 

Figures 6 and 7 are plots of slip coefficient against mean free 
molecular path for both nitrogen and air. The results with nitrogen 
appear to bear out the reciprocal density for the hl^er-temperature 
isotherms, but the 78° Isotherm falls below the others. The 80° iso- 
therm for air also falls below the straight line drawn for the 306 ° 
and 275° Isotherms, and the lowest pressure point at 19^° also falls 
low. 


Region of -Free-Molecule Viscosity 

Theory (see reference 1, p. 300) indicates that, at pressures 
sufficiently low to permit neglect of collisions between molecules , the 
viscous drag exerted by a gas enclosed between two parallel plates 
moving with a relative velocity U is -ZU. The value Z is the free- 
molecule viscosity 6md is represented by 

Z = P/(2jtRT)l/2 (9) 


where R and T represent the molecular gas constant and the absolute 
ten^erature, respectively. This expression indicates that the drag 
should be directly proportional to P and inversely proportional to the 
square root of T. Figure 8 Indicates that this relationship was reached 
at pressures a little below 0.02 millimeter of Eg for both air and 
nitrogen. 


The Ohio State University 

Columbus, Ohio, November 20, 1950 
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TABLE I 


APPARATUS AND WIRE CONSTANTS FOR MANGANIN SUSPENSIONS 


Wire mmi'ber 

Wire constant 

Apparatixs consteint 

1 

0. 00036 

5 . 35 ^^ 

2 

. 00031 

5 . 3 l ^5 

3 

.00021 

5.361 
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TABLE II 

VISCOSITY OF NITROGEN AT VARIOIB TEMPERATURES 


Pressure 

T] 

Pressure 


(inm Hg) 

(poises) 

(mm Hg) 

(poises) 

306.05° K 

273.25° K 

503.0 

1809.2 X 10-T 

550.0 

1661.7 X 10"'^ 

310.0 

1809.2 

403.0 

1661.0 

100.0 

1807.1 

98.0 

1659.4 

68.0 

1804.8 

70.0 

1658. 3 

9.62 

1788. 2 

6.66 

1637.2 

5.07 

1771.5 

3.22 

1610.2 

2.38 

1735.5 

2.34 

1595.4 

1.15 

1651. 3 

1.12 

1515.6 

.574 

1525.7 

.485 

1398.7 

.243 

1308 . 4 

.220 

1249.8 

.100 

922.2 

.097 

937.2 

.056 

6 l 4.2 

.034 

445.4 

.027 

359.4 

.022 

325. T 

.012 

199.3 

.010 

185.6 

.0062 

103.3 

.0030 

53.4 

.0010 

23.4 

. 00041 

9.5 

. 00017 

6.75 

. 000026 

3.1 

193.90° K 

78.60° K 

■ 417-5 

1258. 7 

295.0 

549.5 

102.0 

1257.6 

104.0 

549.5 

72.5 

1256.3 

70.0 

550.2 

9.04 

1250.3 

8.56 

548.1 

4.32 

1243.3 

4.76 

547.2 

1.98 

1227.0 

2.60 

544.2 

.816 

1177.8 

1-17 

541.2 

.353 

1079. 0 

.440 

540.1 



.220 

524.4 



.132 

516.2 
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TABLE III 

VISCOSITIES OF AIR AT VARIOUS TEMPERATURES 


Pressxire 
(mm Eg) 

T 1 

(poises) 

Pressure 
(mm Eg) 

T1 

(poises) 

306. 05° K 

275.00° K 

696. 0 

1880, 9 X 10"^ 

195.0 

1735.4 X 10"^ 

402.0 

1879. 4 

98.0 • 

1735.2 

100.0 

1877.4 

50.5 

1724. 3 

70.0 

1876. 0 

26.5 

1724. 9 

9.70 

1857 . 6 

8.38 

1712.8 

5.38 

1838 . 3 

4.46 

1698. 0 

1.76 

1751. 1 

.880 

1533.0 

.440 

1459. 3 

.507 

1351. 8 

.220 

1107. 9 

.232 

1046.9 

.088 

802.9 

.088 

770.0 

.035 

458.4 

.044 

537.0 

.020 

300.0 

.022 

338.4 

.010 

161.9 

.0093 

167.7 

.0042 

66.9 

.0042 

73.7 

.00054 

5.7 

. 00033 

10.0 

193.90° K 

80.30° K 

495.5 

1309.2 

314.0 

575.7 

105.0 

1304. 1 

205.0 

579.2 

70.0 

1303.8 

49.0 

581.6 

7.72 

1290,6 

25.0 

581.7 

3.57 

1279.4 

9.35 

582.2 

1.76 

1263.4 

4.68 

579.9 

.970 

1239.2 1 

2.18 

579.8 

' .308 

1150, 2 

• 1.00 

571.4 


f 

.397 

547.0 



.243 

542.5 


















12 


MCA TN 2546 


TABLE IV 

SLIP COEPFICIEIJTS OF NITROGEN AT VARIOUS TEMPERATURES 
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TABLE V 

SLIP COEFPTCIEMTS OF AIR AT VARIOUS TEMPERATURES 
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TABLE VI 

VALUES OF CONSTANT k IN EQUATION (8) 


Temperatiire 

Constant k for - 

(°K) 

Nitrogen 

Air 

306.05 

4 . 81 X lo 2 

5.90 X 102 

275.00 

273.25 

193.90 

3.96 

2.92 

5.50 

2.37 

80.30 

78.60 

.53 

1.08 
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Figure 1.- Viscosity of nitrogen against logarithm of pressvire at 

fixed teraperatrares . 
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Figure 2.- Viscosity of air against logarithm of pressure at 

fixed temperatures. 
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Figure 3-- Slip coefficient of nitrogen to 0.1 millimeter of mercury 
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Figure 4 .- Slip coefficient of air to 0.1 millimeter of mercury. 
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Figiore 6.- Slip coefficient of nitrogen against mean free path. 
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